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We study the conditions for successful Affleck-Dine baryogenesis and the origin of gravitino dark matter in 
GMSB models. AD baryogenesis in GMSB models is ruled out by neutron star stability unless Q-balls are 
unstable and decay before nucleosynthesis. Unstable Q-balls can form if the messenger mass scale is larger than 
the flat-direction field <J> when the condensate fragments. We provide an example based on AD baryogenesis 
along a = 6 flat direction for the case where mj,i2~2 GeV, as predicted by gravitino dark matter from Q-ball 
decay. Using a phenomenological GMSB potential which models the <J> dependence of the SUSY breaking 
terms, we numerically solve for the evolution of <I> and show that the messenger mass can be sufficiently close 
to the flat-direction field when the condensate fragments. We compute the corresponding reheating temperature 
and the baryonic charge of the condensate fragments and show that the charge is large enough to produce late- 
decaying Q-balls which can be the origin of gravitino dark matter. 
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I. INTRODUCTION 

Supersyrmnetry (SUSY) is widely viewed as the most likely 
theory beyond the Standard Model (SM). In SUSY exten- 
sions of the SM, the spectrum of sparticle masses is strongly 
dependent on the mechanism for SUSY breaking. Gauge- 
and gravity-mediated SUSY breaking are the most commonly 
considered mechanisms. In the case of gauge-mediated SUSY 
breaking (GMSB), the gravitino can be much lighter than the 
mass of the sparticles, making it a natural candidate for the 
lightest SUSY partner (LSP) and dark matter For a complete 
cosmology of GMSB, we also need a mechanism for baryo- 
genesis. In the MSSM, Affleck-Dine (AD) baryogenesis 1 1] is 
a particularly simple and effective way to generate the baryon 
asymmetry. A complex flat direction field gains a large ex- 
pectation value, which later begins coherent oscillations in the 
real and imaginary directions. The interaction of the evolving 
field with B violating operators induces a net asymmetry in 
the flat-direction condensate. 

It is well-known that AD condensates are unstable with re- 
spect to spatial perturbations |2-9]. These grow and fragment 
the condensate, producing Q-balls. The cosmology of Q-balls 
depends on the dimension of the B-violating operator and the 
nature of SUSY breaking. In gravity-mediated models the Q- 
balls are unstable and, if their charge is large, these may de- 
cay below the freeze-out temperature of neutralino dark matter 
llHH] - In this case dark matter and baryon number are simul- 
taneously produced when the Q-balls decay, with yidm = 3/1^ 
by R-parity conservation. Therefore Q.b/^dm = iih,/ {3moM), 
where m„ is the nucleon mass. For neutralino dark matter 
with mass ^100 GeV, as expected in gravity-mediated SUSY 
breaking models, this produces a value much smaller than the 
observed baryon-to-dark matter ratio w 1 /6, which requires 
moM ~ 2 GeV. In this case the neutralino DM candidates 
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must annihilate rapidly after production, pointing to Higgsino- 
like dark matter [10]. 

While this is an appealing mechanism for non-thermal dark 
matter, the possibility of directly producing both dark matter 
and the baryon asymmetry from Q-ball decay is lost. This is 
only possible with a dark matter candidate of mass « 2 GeV. 
In [11] a light axino LSP in gravity-mediated SUSY breaking 
was proposed. Alternatively, in GMSB models a natural low- 
mass DM candidate is the gravitino. 

However, Q-ball formation in GMSB AD baryogenesis is 
difficult to achieve without phenomenological problems lfT2l - 
I14fl . The flat-direction potential depends on the messenger 
mass. For field values much larger than the messenger mass, 
the potential becomes almost flat due to the suppression of 
gauge-mediation by the flat-direction field. If the magni- 
tude of the flat-direction field when Q-balls form is much 
larger than the messenger mass, then the Q-balls will have an 
energy-per-charge, /i/Q, which decreases as Q^^^'^ |2]. For 
large enough Q, E /Q is less than the nucleon mass, imply- 
ing that the Q-balls are absolutely stable if they carry baryon 
number This leads to severe problems. A stable Q-ball cap- 
tured by a neutron star will absorb baryon number, reducing 
its mass and eventually destabilizing the neutron star This 
severely constrains AD baryogenesis in such GMSB mod- 
els, ruling out all flat directions lifted by conventional B- 
conserving non-renormalizable potential terms (l^. More 
precisely, fi/g °^ fi"'/'* is true if the AD field (|) in the Q-ball 
is on the approximately constant part of the potential. (These 
are referred to as flat direction (FD) Q-balls in [12].) For suf- 
ficiently large Q-ball charge, (|) becomes large enough for the 
non-renormalizable terms to become important in the Q-ball 
solution. This happens once Q >Qc, where for the d — 6 
Q-balls we consider Qc ~ lO'^^ 01211 . For larger charge, the 
Q-balls become thin-walled, with constant field magnitude 
(^c ^ lO'"* GeV inside and radius increasing as Q increases. 
(These are referred to as curved direction (CD) Q-balls in 
il2ll .) If Q-balls have a density comparable to dark matter (as 
we would expect if the Q-balls are stable and the baryon asym- 
metry is due to Affleck-Dine baryogenesis), then they will be 
excluded by neutron star stability if the non-renormalizable 
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lifting terms conserve baryon number. FD type Q-balls im- 
ply a neutron star lifetime ~ 10^" years, whereas CD Q-balls 
imply a lifetime ~ 1500 years lll2ll . The initial charge of Q- 
balls from AD baryogenesis will typically be less than Qc and 
therefore they will initially be FD type Q-balls. However, 
since Qc <C Qns, where Q„s ^ 10^^ is the baryon number of 
a neutron star, the charge of any FD Q-ball trapped by a neu- 
tron star will rapidly grow by baryon absorption to become 
larger than Qc (for a constant absorption rate, this will occur 
in a time ^ Qc/Qns x 10^" years ^ 10^^^ s), transforming 
them in CD Q-balls and destabilizing the neutron star. 

AD baryogenesis in GMSB could work if the messenger 
mass is sufficiently large that condensate fragmentation oc- 
curs when the amplitude of the flat-direction scalar is in the 
approximately quadratic part of the potential. In this case 
the resulting Q-balls will be unstable. It may then be possi- 
ble to produce dark matter via late decay of Q-balls to out- 
of-equilibrium NLSPs which subsequently decay to gravitino 
with mass m^i2 ~ 2 GeV.' 

In the case of gravity-mediated SUSY breaking, it is known 
that d = () flat directions lifted by superpotential terms of the 
form [u^d'^d'^)^ can produce Q-balls with large enough baryon 
number to decay after neutralino freeze-out but before nucle- 
osynthesis ||3j, i4j]. We may expect similar behaviour in the 
case of GMSB flat directions with sufficiently large messen- 
ger mass. Therefore in the following we will focus on the 
example of d = 6 flat directions ^. Our method will be to nu- 
merically study a phenomenological GMSB flat-direction po- 
tential which models the suppression of the SUSY-breaking 
masses at field values greater than the messenger scale. 

Our paper is organized as follows. In Section 2 we intro- 
duce the phenomenological GMSB flat-direction potential. In 
Section 3 we consider AD baryogenesis and condensate frag- 
mentation along a c/ = 6 MSSM flat-direction in a GMSB 
model with m^ji ~ 2 GeV. We solve the field equations to 
obtain the baryon density and so determine the reheat temper- 
ature. We then perform a semi-analytical study of condensate 
fragmentation using the phenomenological potential and esti- 
mate the time of formation and baryonic charge of the con- 
densate fragments. In Section 4 we present our conclusions. 



II. FLAT-DIRECTION POTENTIAL IN GMSB 

GMSB models are based on SUSY breaking in a hidden 
sector which is transmitted to the MSSM via vector pairs of 
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messenger fields with SM gauge charges. The messenger 
superfield scalar components acquire SUSY breaking mass 
splittings from their interaction with the hidden sector. The 
messengers then induce masses for MSSM gauginos at 1-loop 
and soft SUSY breaking scalar mass squared terms at 2-loops. 

This is true so long as the flat-direction field <I> does not 
give the masses to the gauge fields which are greater than 
the messenger mass M,,,. Once |<I>| ^ M,„, the transmission 
of SUSY breaking via gauge fields is suppressed^. In gen- 
eral, all SUSY breaking mass terms and A-terms are expected 
to be proportional to the order parameter of SUSY breaking, 
< Fs>, therefore once I'Pl^ M„, all SUSY breaking mass pa- 
rameters are expected to be proportional to <Fs > /\^\ lUTIl . 
Thus we must include a lO]^^ suppression factor for the soft 
SUSY breaking mass terms at large |<I>|. A 2-loop calculation 
shows that the flat-direction effective potential has in fact a 
squared logarithmic dependence on |<I>| iHtIi . 

Therefore to model d = 6 AD baryogenesis in GMSB 
models, we consider the following phenomenological flat- 
direction potential 



M,„ 



1 +K\n 



1^ 
Mi 



'3/2 



1 +K\n 



M2 



where 



{AW + h.c: 



w = 



(1) 



(2) 



We include the factor 6! in Eq. (|2]) so that the physical strength 
of the interactions is dimensionally of the order of M. The 
first term in Eq. ^ is due to GMSB with messenger mass 
M,„. The factor multiplying this takes into account 1-loop ra- 
diative corrections due to gaugino loops once |<I>| "z M,„, with 
Kk, —(0.01 —0.1) [3-5]. The second term is due to gravity- 
medated SUSY breaking including the 1-loop correction term 
K. (For simplicity we set K ^ K.) We have also included a 
Hubble correction to the mass squared term. We do not in- 
clude a Hubble correction to the A-term, as the interaction 
leading to such terms is typically proportional to the inflaton 
field and so averages to zero over its coherent oscillations list 
The scale of the non-renormalizable terms, M, will depend on 
the origin of the superpotential term responsible for lifting the 
flat direction. This could be gravitational, M ^ Mp (where 
= 2.4 X 10'** GeV) or due to the new physics suggested by 



' More accurately, g\^\ ^ M,,,, where g is the relevant gauge coupling. For 
the (u'^d'^d'^)^ direction of interest to us here, g is the strong coupling. Since 
g~ I, for simplicity we will not include g explicitly in the potential. 
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MSSM coupling constant unification, M ^ Mu w 10 ° GeV. 
In the following we will focus on the case where M ^ Mp. For 
the A-term we consider 



'«3/2 - 



1/2 • 



(3) 



Therefore the messenger mass is related to mj and by 

^2 V3K7W3/2Mp 



Thus 



1671:2 



(10) 



The first term in Eq. (|3]l represents the A-term due to gravity- 
mediated SUSY breaking. The second term models the A- 
term in gauge-mediated SUSY breaking at |<I>| M,„, which is 
generated at 1-loop from the gaugino masses and is therefore 
suppressed relative to the A-term in gravity-mediated models, 
with fl„ ^ 0.01. The suppression factor (1 + |<I>|2/m2)-i/2 
models the 1 /|<I>| suppression of the GMSB A-term at |<I>| > 

As a reference model for GMSB we will consider a 
metastable SUSY breaking sector lfl9l] . although our results 
depend only on the SUSY breaking contribution to the mes- 
senger scalar masses and therefore can easily be adapted to 
GMSB models more generally. The superpotential of the 
SUSY breaking sector is 



(4) 



where S is the SUSY breaking field and /, / are the messenger 
fields. The Kahler potential is 



A4 



(5) 



which is the generic form of Kahler potential for GMSB in a 
metastable vacuum. The potential is then 



y = |A<2-K//|2(i + n^ + o(^ 



\sl 

A2 



\s\' 

A4 



|k5/ + M/|2 + |k5/ + M/|2 



(6) 



There is a local SUSY breaking minimum at 5 = / = / = if 
m2 > K^2^ jjj which case Fs = fr. At the local SUSY break- 
ing minimum the SUSY messenger mass and SUSY breaking 
mass squared splitting are 



M,„ = M ; F,„ = kFs = k/j- 



(7) 



The masses squared of the scalar components of the mes 
sengers /, / are therefore m2 — nij = m2 ± yqr-, while ths 
fermion mass is M,„. Soft SUSY-breaking masses for gaugi 
nos and scalars in the MSSM sector of order are then gen 
erated by messenger loops 



167l2 M„ 



167t2 M,n 



(8) 
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From this we see that a large messenger mass scale is required 
to have m3/2 ~ 2 GeV. The importance of this for d — 6 AD 
baryogenesis is that, for plausible values of M, the messenger 
mass can be comparable to the initial value of the flat direction 
scalar when <I> begins to oscillate and the baryon asymmetry 
is generated. To get a rough estimate of the value of |<I>| at 
the onset of oscillations, lOlojc, we compute the value of the 
minimum of y(<I>) for the case where only the — c//2|<J)|2 term 
and the SUSY non-renormalizable term are included. Setting 
c//2 = m2, the value when the potential is destablized by the 



mass squared term, then gives 



If la«. ^ -T- 



5!2 
5 



3x lO^^GeVx 



V 100 GeV/ V 



1018 GeV J 



(12) 



The value of M is a free parameter in AD baryogenesis mod- 
els. A natural possibility is that the non-renormalizable terms 
are associated with the completion of the theory at the Planck 
scale. However, even with this assumption there are differ- 
ent possibilities for how the Planck scale enters. If one ex- 
pects the strength of the interaction to be dimensionally de- 
termined by Mp, then the factorial term should be included 
in the superpotential and so M sa Mp. But it is also of- 
ten assumed that the Planck scale enters in the Lagrangian 
rather than the vertex, in which case W = <I>Vm^ and so 

M w Mp/(6!)i/3 = 2.7 X 10'^ GeV. We will therefore con- 
sider a range of M from 2 x lO'^ GeV to Mp'^. 

For the case where ^2^^ 1, 100 GeV and M = 

2 X lO'^ GeV, the messenger mass scale is 5 x lO'-' GeV 
while \^\osc ~ 8 X 10^^ GeV. Therefore it is possible that the 
messenger mass scale can be close to the value of |<1>| at the 
onset of oscillations. In addition, there will be a significant 
time delay between the onset of oscillations and the eventual 
fragmentation of the AD condensate and Q-ball formation. 
Therefore even if |<I>| is greater than M,„ when oscillations 
begin, it is possible that |<I>| will be less than M„, when the 
condensate fragments. If Q-balls form when |<I>| "z M,„, then 



The gravitino mass is 



'«3/2 



Fs 



V3Mp V3Mp 



Such scales may also arise in the case where M is due to exchange of 
heavy particles of mass My ~ lO'* GeV, once suppression of the non- 
renormalizable operators due to couplings is included. 
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E /Q will not be strongly suppressed, in contrast to the case of 
conventional GMSB Q-balls with |4>| > M„,. 

Thus d = 6 flat directions, combined with /2 ~ 2 GeV 
andM« 10^^ - lO'** GeV, may lead to condensate fragmen- 
tation in a region of the potential where unstable Q-balls can 
form. In order to better understand the evolution of the flat- 
direction condensate, in the next section we will numerically 
evolve <I> using the phenomenological GMSB potential given 
in Eq. ([T]). 



III. AFFLECK-DINE BARYOGENESIS AND 
CONDENSATE FRAGMENTATION 

A. d = 6 AD Baryogenesis in GMSB with a large messenger 
mass 
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FIG. 1 : Growth of the baryon asymmetry per comoving volume as a 
function of |<I>|/M,„. 



The baryon asymmetry in the late-time coherent oscilla- 
tions of the real and imaginary parts of <I> is induced by the 
B-violating A-term. The baryon number density is ng = ng/3, 
where ng is the global t/(l) charge density corresponding to 
the number asymmetry in <I> particles: 



Q ^ / ( - 4)t(j) j = (|)j(j)2 _ (i)i(|)2 



(13) 



The factor 1/3 accounts for the baryon number of <I> along the 
{u'^d'^d^')^ direction. In the following we will use subscript b 
to denote quantities at the time when the baryon asymmetry 
in a comoving volume becomes fixed. The present baryon 
asymmetry to entropy is then 



(14) 



The observed baryon asymmetry is {ns/s)oh.<i ~ 1.8 ±0.1 x 
10^^". This allows us to fix the reheating temperature for a 
given flat-direction potential by computing Hi, and nsh- 

In Fig.l we show the growth of the baryon asymmetry 
per comoving volume as a function of |<I>|/Mm for the case 
M,„ = 5 X 10'^ GeV, m, = 100 GeV, K = -0.1 and M = Mp, 
where |<t>| is the maximum value of the magnitude of <I> on the 
ellipitical trajectory. (We assume a„ = 0.01 in all our calcula- 
tions.) This shows the rapid increase in the asymmetry once 
the <I> oscillations begin, with the asymmetry becoming fixed 
once |<I>|/M,„ 3. In Fig. 2 we show the trajectory of <I> in the 
complex plane when |<I>|/M„, = 0.54. A noteworthy feature 
of this is that the trajectory is not yet ellipitical but in fact pre- 
cesses, even though the baryon asymmetry per comoving vol- 
ume is constant. This indicates that the GMSB deviation from 
a jOp potential is still significant at this time, even though the 
A-terms have become negligible and the baryon number in a 
comoving volume is constant. 

In Table 1 we give the AD baryogenesis parameters and Tr 
for a range of M and M„j for the case = 100 GeV. In Table 
2 we show the case with nis = 500 GeV. From this we see that 
in most cases the baryon asymmetry forms when |<I>|/M,„ is 
larger than or of the order of 1 . 
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FIG. 2: Trajectory of <I> in the complex plane after formation of the 
baryon asymmetry. The ellipse precesses but the baryon asymmetry 
per comoving voume is constant. 



B. Perturbation growth and condensate fragmentation 



To study the growth of perturbations of the AD condensate 
we will use the approach of |2]. This assumes that the baryon 
asymmetry in the condensate is maximal i.e. the trajectory of 
<I> is circular in the absence of expansion ("Q-matter"). This is 
not true for GMSB models with a„ ^ 1 and /2 ^ which 
produce ellipitical trajectories. However, the growth rate of 
energy density perturbations for the ellipitical trajectory will 
be similar, since in this case the growth of perturbations will 
be determined by the negative pressure of the condensate av- 
eraged over oscillation cycles, which turns out to be essen- 
tially independent of the field trajectory for a given potential 
and maximum value of the field on the trajectory, (^i„ax- The 
average pressure for the limiting case of a linear oscillating 
condensate (i.e. a real field) has been discussed in 1201, where 
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M 


Mm 


nsb 


Hh 


if) 

M„, 


Tr 


100 


2.4 X 10"^ 


5 X 10'^ 


1.4 X 10^*^ 


6.8 


2.5 


14.4 


100 


2.4 X 10'^ 


5 X lO'^ 


4.3 X 10^^ 


1.6 


22.6 


2.4 


100 


2.4 X 10"^ 


1 X 10'2 


4.3 X 10^5 


0.4 


75.0 


14.7 


100 


2.0 X 10'^ 


5 X 10'^ 


3.7 X 10^'' 


15.4 


0.5 


2701 


100 


2.0 X 10'^ 


5 X 10'2 


2.9 X 10^^ 


5.7 


14.1 


464 


100 


2.0 X 10'^ 


1 X 10'2 


1.2 X lO^* 


2.0 


17.3 


140 



TABLE I: AD baryogenesis parameters for m, = 100 GeV. (Dimensionful 
quantities in GeV units.) 





M 


M,„ 


nBh 


Ht 


1*1 


Tr 


500 


2.4 X 10"* 


1 X 10'^ 


1.2 X lO^s 


9.7 


14.9 


33.3 


500 


2.4 X 10"* 


1 X 10'2 


6.7 X 10^^ 


1.9 


119.6 


2.1 


500 


2.0 X 10'^ 


1 X 10'^ 


9.3 X 10^^ 


33.0 


3.2 


5.1 X 10^ 


500 


2.0 X 10'^ 


1 X 10'2 


1.2 X 10^^ 


6.8 


19.7 


845 



TABLE II: AD baryogenesis parameters for nij = 500 GeV. (Dimensionful 
quantities in GeV units.) 



the average pressure was shown to be 



P^ + l=2^^ ^'""^ ■ (15) 



< 



Perturbations of the energy density of wavenumber k then 
grow according to 



8pk 



<p> IkpSpk 



(16) 



In this we have neglected expansion since the growth rate is 
typically large compared with H. In the case of the circular 
condensate, the pressure is constant and determined by the 
magnitude of the field and the potential, 



2 



P K. 

2 



(17) 



V 



This expression allows us to derive simple expressions for the 
pressure and growth rate of perturbations which depend only 
on the potential and its derivatives, unlike the linear oscilla- 
tion case where the average pressure must be calculated nu- 
merically via Eq. ( fTSl l. The key observation that allows us to 
also use Eq. ST% for non-circular trajectories is that the pres- 
sures in a circular condensate and in a linear condensate for 
a given (^,„ax are essentially identical. For the case of a poly- 
nomial potential, the pressures can be calculated exactly in 
both cases and are exactly equal, < p > /p = (n — 2)/(n + 2) 
for y('t>) o<: |<I>|" ||20^. For a non-polynomial potential there 
is no analytical expression and the integrals in Eq. ( fTSl ) must 
be calculated numerically. (This is not straightforward as 
the integrand in the denominator of Eq. ( fTSl ) is divergent as 



V Vmax-) We have numerically integrated Eq. ( fTsT i for the 
potential Eq. ([T]| and find that the negative pressure is equal 
to within 10% for a given (|)„,„r. Therefore the growth rate 
of perturbations of the circular and linear condensate will be 
essentially the same. We can expect the same to be true for 
ellipitical trajectories which are between these limiting cases. 
This allows us to use the more convenient expression for the 
pressure Eq. ( [TtT i to calculate the growth of perturbations of 
the ellipitical condensate and so the time when they become 
non-linear The similarity of the growth of perturbations in 
the circular and linear condensates is physically reasonable, as 
the growth of perturbations is essentially due to an attractive 
interaction between the scalar particles, and the strength of 
the interaction is the same for real and complex scalars which 
have the same potential. This approach is also in agreement 
with the numerical results of [8], where it was found that the 
size of the condensate fragments in the circular and elliptical 
cases differs by a factor of at most 2. 

We first consider condensate fragmentation for a circular 
condensate in the potential Eq. ([T]). The flat-direction field 
can be written as 



4>= ^7;(x,f)e'"M 
v2 



The <I> field equations are then 



1 



R 



Q. + 3HQ.- + 

R 



j^R 



VQ..VR ^ 



(18) 



(19) 



R + 3HR 



rV^R 



Q?R 



;{VafR + V =0, (20) 



where V = dV /dR. In this it has been assumed that V is a 
function of R alone, so these equations can be used to study 
growth of perturbations only once the A-terms are no longer 
significant in V i.e. once H < Hi,. The growing mode solu- 
tions have the form 



(21) 



Substituting these into Eq. ( fTsT i and Eq. iT% and assuming 
that a = S = constant gives ||2(] 



a 



fl2 



2R 

3Ha + ^ + —a 



R 



a" 



3//aH 



'Q.'^ + U {R) 





R 











a = , 



(22) 



a is assumed to be constant because the rate of growth of 
perturbations is large compared with the rate of expansion of 
the Universe, in which case the time dependence in Eq. ( l22b . 
which is entirely due to expansion in the case of circular tra- 
jectory, can be neglected. In this case R can also be set to zero, 
as this is non-zero only because of expansion. In this case the 
equation simplifies to 



a 



fl2 



a 



-^-n^ + y (R) 



-4aV=0. (23) 
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This can be solved for a 



a" 



(y" + 3a2) 



Q.^ - V" 



16- 



«^ (y" + 3n2)^ 



(24) 

From this the largest value of |k| for which growth is possible 
is 



V"+3Q? 



(25) 



where we have defined Km — \kmax/a\. For a homogeneous 
Q-matter background with R ^ R = 0, Eq. ( l20l i implies that 

The solution for a in Eq. ( l24b gives the value of a(f ) on a 
time-scale short compared to H^. One can then re-introduce 
the t dependence of R and Cl due to the expansion of the Uni- 
verse. The growth factor ^(f ) is then given by 



Sit) 



a{t)dt 



(26) 



where is the time at which a mode of wavenumber k begin 
to grow. 

For a given mode k, the perturbation of the condensate will 
only begin to grow once |k/fl| < K,„. Suppose for a given 
mode this occurs at a time . The subsequent growth of the 
perturbation follows from 



Ja, aH 



(27) 



where 



a(k,a) 



fl2 (y" + 3a2) 



1/2 



(28) 



With H °^ a ^/2 during inflaton-domination, we find from 

Eq. dan 



5(k,fl) = 2a(k,fl) (l - ^\ H 



(29) 



At a given value of a, the mode with the maximum growth is 
found by maximizing 5(k, a) with respect to |k|. Since a o= |k| 
and a^,{k) = A',^' |k| we find 



5(k) oc |k| 



ML 



(30) 



This is maximized at 



a 9 



(31) 



Therefore at each value of a, the mode which has the largest 
growth will have wavenumber given by Eq. ( [3TT l. The value of 
S for this mode is 



(32) 



The condition for the condensate to fragment is 

1 , (33) 



dR _ 5Ro s > 



R Ro 

for the mode with wavenumber Eq. ( l3Tl l. The diameter of the 
fragments is then given by 



"kfn 



2n 



{\^\frag/a) 



(34) 



while their global f/(l) charge (Q — 3B) is 

3 



471 / " kfrag \ ' 



(35) 



We expect the same value of S to apply to the elliptical tra- 
jectory when R = Rmax, the maximum value of R for the el- 
lipse. We will therefore use this to estimate the time of frag- 
mentation and the value of |<I>|/M„, at fragmentation. Given 
the energy density and baryon density when the condensate 
fragments, we can then compute the energy and baryon num- 
ber of the fragments. 

To compute the time of fragmentation of the condensate, 
we need a spectrum of initial perturbations. These are related 
to the primordial density perturbations. For the initial pertur- 
bation entering the horizon we expect dR/R^ 5i2 /i2 ss 5p/p, 
the usual primordial density perturbation. This is certainly 
true of the amplitude, since R oc at the minimum of the 
potential due to the Hubble mass correction. For the phase, 
once the perturbation enters the horizon we can expect field 
dynamics to transfer the perturbation of R into a correlated £2 
perturbation of a similar magnitude. There could also be addi- 
tional phase fluctuations independent of the amplitude which 
will result in isocurvature baryon perturbations, but these are 
model dependent and we will not consider them here. With 
5Ro/Ro ~ 5p/p ~ 10^^, the condition for fragmentation be- 
comes 5 11. 

In Figure 3 we show the growth factor 5 as a function 
of \<P\/M,„ for the case m, = 100 GeV, K = -0.1, M,„ = 
5 X lO'"* GeV and M = Mp. It can be seen that fragmenta- 
tion occurs at |<I>|/M,„ « 0.1, even though the baryon asym- 
metry forms at a much larger value, |4>|/M,„ « 2.5. In Ta- 
ble 3 we show the condensate fragmentation parameters for 
the case nis — 100 GeV and in Table 4 for the case where 
nis = 500 GeV. We see that fragmentation can easily occur 
with |4>|/M„, < 1, in particular if M < Mp and \K\ < 0.1. 

The field trajectory at fragmentation is strongly ellipitical, 
indicating that the total energy density in the <I> field is large 
compared with the energy density in Q-matter of the same 
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charge density. We can estimate this ratio for a strongly ellip- 
tical trajectory by 



p vm) 

{E/Q)nQ m^riQ 



(36) 



This is correct for the case where |<I>| for the Q-matter config- 
uration is sufficiently less than M,„ that it can be considered 
to be described by |<I>p potential with E/Q^i m^. Otherwise 
E /Q will be less than nis, so that Eq. ( l36b will be an underesti- 
mate of rE- rE also gives the ratio of the energy of the conden- 
sate fragments to the energy of a Q-ball of the same charge. 
For the cases being considered we find that rs is mostly in the 
range 10-60 for m, = 100 GeV and 20-90 for m, = 500 GeV. 

We next estimate the decay temperature of the Q-balls 
under the assumption that the condensate fragments evolve 
into gravity mediated-type Q-balls of the same charge, corre- 
sponding to Q-balls in a potential V = m^\<P\^{l +Kln{\<P^). 
The decay rate of a Q-ball has an upper bound given by ll2lll 



(37) 



where the Q-ball solution is proportional to exp{mt) and A 
is the area of the Q-ball. The upper limit is expected to be 
saturated for flat-direction Q-balls. The decay temperature is 
then 



1/2 



(38) 



where kj = {TZ^g{T)/90y/^ « 4.7 using g{T) « 200 for the 
MSSM. Dimensionally we expect CO ^ m and R ^ so 
we can parameterize (sy'R^ as (a^R^ = fgin. In particular, for 
gravity-mediated Q-balls co w m and R^ = 2/ {\K\m^), there- 
fore fg = 2/\K\. This assumes <I> decay to scalar pairs is 
kinematically excluded, otherwise the decay rate can be en- 
hanced by a factor /j ~ 10^ since the decay is not Pauli- 
blocked within the volume 10] • The decay temperature can 
then be written as 



Td = 60 MeV/g' 



1/2 



100 GeV 



1/2/10- 



,20 \ 1/2 



Q 



(39) 



In Table 5 we show the gravity-mediated type Q-ball prop- 
erties for the case m, = 100 GeV and in Table 6 for the case 
nis — 500 GeV. In most cases the Q-ball decay temperature 
is well above the temperature of the onset of nucleosynthe- 
sis « 1 MeV but well below the typical temperature of NLSP 
freeze-out w Mesp/2Q ^ 5 GeV. (The decay temperature will 
be somewhat higher if fs > 1 .) Therefore such Q-balls can be 
a source of non-thermal NLSPs which subsequently decay to 
gravitino dark matter 

The condensate lumps will evolve into gravity-mediated Q- 
balls if at the end of their evolution |<I>|/M,„ <C 1, so that the 
dynamics of the Q-ball are entirely due to the radiatively- 
corrected |<I>p potential. The initial values of |<I>|/M„, from 
Tables 3 and 4 are not always very much less than 1 . However, 



since the energy in the initial condensate fragments is larger 
than the final energy in the Q-ball by a factor r^, the field 
in the final Q-ball configuration will be smaller by a factor 
~ 1 1 \frE- In most cases this is enough to ensure |<I>|/M,„ < 1 
for the Q-bafls. 

Note that if |<I>|/M,„ is less than 1 but not very small, the 
Q-ball may correspond to a new type of Q-ball which is inter- 
mediate between the gravity- and gauge-mediated cases. We 
will study the structure of such intermediate Q-balls in a fu- 
ture discussion [i22il . 



rris 


K 


M 


M,„ 


1*1 




"Bfrag 


rE 


100 


-0.1 


2.4 X 10'^ 


5 X 10'^ 


0.12 


54.4 


5.0 X 10^^ 


19.1 


100 


-0.01 


2.4 X 10'^ 


5 X 10'^ 


0.017 


21.1 


2.9 X lO^"* 


24.3 


100 


-0.1 


2.4 X 10'^ 


5 X 10'2 


3.2 


24.9 


3.0 X 10^^ 


5.8 


100 


-0.1 


2.4 X 10'^ 


1 X 10'2 


12.3 


9.6 


4.3 X 10^^ 


1.6 


100 


-0.1 


2.0 X lO''^ 


5 X 10'^ 


0.006 


45.1 


1.8 X 10^^ 


43.3 


100 


-0.01 


2.0 X lO''^ 


5 X 10'^ 


0.0009 


14.6 


1.1 X 10^' 


57.9 


100 


-0.1 


2.0 X 10^'' 


5 X 10'2 


0.26 


55.1 


2.7 X lO^^^ 


16.3 


100 


-0.1 


2.0 X lO''^ 


1 X 10'2 


2.3 


29.5 


7.0 X lO^^^ 


7.3 



TABLE III: Condensate fragmentation parameters for = 100 GeV. (Di- 
mensionful quantities in GeV units.) 





K 


M 


Mm 


1*1 


Km 


"Bfrag 


rE 


500 


-0.1 


2.4 X 10'^ 


1 X 10'^ 


2.2 


149.9 


7.3 X lO^'' 


31.8 


500 


-0.01 


2.4 X lO'** 


1 X lO'^ 


2.2 


153.3 


7.3 X 10^^ 


31.8 


500 


-0.1 


2.4 X 10"^ 


1 X 10'2 


21.2 


30.6 


8.3 X 10^^ 


1.9 


500 


-0.1 


2.0 X 10'^ 


1 X 10'^ 


0.15 


276 


4.7 X lO^"* 


72.2 


500 


-0.01 


2.0 X 10'^ 


1 X 10'^ 


0.03 


118 


1,2 X 10^^^ 


85.3 


500 


-0.1 


2.0 X 10'^ 


1 X 10'2 


3.8 


111 


2.1 X 10^5 


20.0 



TABLE IV: Condensate fragmentation parameters for ; 
mensionful quantities in GeV units.) 



: 500 GeV. (Di- 



0.2 0.4 0.6 0.8 

3)/ Mm 



FIG. 3: Growth factor S(0 as a function of |<5|/Af„; 
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K 


M 


M,n 


Q 


rrf(MeV) 


100 


-0.1 


2.4 X 10'** 


5 X 10'^^ 


1.1 X 10^^^ 


8.1 


100 


-0.01 


2.4 X 10'** 


5 X 10'^^ 


4.0 X 10^2 


42 


100 


-0.1 


2.4 X 10'** 


5 X 10i2 


8.0 X 10^4 


1.0 


100 


-0.1 


2.4 X 10'** 


1 X 10'2 


6.5 X 10^5 


0.3 


100 


-0.1 


2.0 X 10'^ 


5 X 10'^^ 


2.6 X lO^" 


167 


100 


-0.01 


2.0 X 10'^ 


5 X 10'^^ 


1.3 X lO^O 


236 


100 


-0.1 


2.0 X 10'^ 


5 X 10'2 


6.5 X 10^1 


33 


100 


-0.1 


2.0 X 10'^ 


1 X 10'2 


9.2 X 10^5 


0.3 



TABLE V: Q-ball charge and decay temperature for = 100 GeV. (Di- 
mensionful quantities in GeV units except Tj.) 





K 


M 


M,n 


Q 


rrf(MeV) 


500 


-0.1 


2.4 X 10'** 


1 X lOi"* 


7.8 X 10^2 


22 


500 


-0.01 


2.4 X 10'** 


1 X 10'^^ 


7.8 X 10^2 


69 


500 


-0.1 


2.4 X 10'** 


1 X 10i2 


1.1 X 10^5 


1.8 


500 


-0.1 


2.0 X 10'^ 


1 X 10'^^ 


7.7 X lO'** 


683 


500 


-0.01 


2.0 X 10'^ 


1 X 10'^^ 


2.7 X 10''* 


3651 


500 


-0.1 


2.0 X 10" 


1 X 10^2 


5.7 X lO^' 


79 



TABLE VI: Q-ball charge and decay temperature for m, = 500 GeV. (Di- 
mensionful quantities in GeV units except Tj.) 



IV. CONCLUSIONS 

In this paper we have revisited Affleck-Dine baryogenesis 
and Q-ball formation in GMSB models. We have been con- 
cerned with two questions: (i) whether unstable Q-balls can 
naturally form in GMSB AD baryogenesis models and (ii) 
whether such unstable Q-balls can be a source of gravitino 
dark matter To answer these questions we have focused on 
a.d ~6 flat direction with m3/2 = 2 GeV, which is consistent 
with gravitino dark matter from Q-ball decay. 

We find that it is possible for condensate fragmentation in 
d = 6 AD baryogenesis to occur when |«1>|/M,„ is less than 
or of the order of 1. The importance of this is that frag- 
mentation does not occur on the plateau of the GMSB po- 
tential and so does not lead to formation of stable Q-balls 
which would destabilize neutron stars. Instead the fragments 
form either in the intermediate region between the \<P\^ region 
and the GMSB plateau or well within the \<P\^ dominated re- 
gion, where we can expect the Q-balls to behave like unstable 
gravity-mediated Q-balls. If the condensate fragments subse- 
quently evolve into Q-balls of the same charge, then the Q- 
balls should be of approximately gravity-mediated form and 
unstable, with decay temperatures typically in the range 10 
MeV - 1 GeV. Decay of such Q-balls to NLSPs could then 
be a non-thermal source of gravitino dark matter. Even if 
I'l'l/^m ^ 1 when the condensate fragments, the final value 
of |<I>|/M,„ in the Q-ball will typically be significantly less 
than 1 once the fragment has lost its excess energy, since the 
energy of the initial condensate fragment is much larger than 
the energy in its baryonic charge. 



In some cases the value of |<I>| will be less than the mes- 
senger scale but not sufficiently small that the Q-ball will be 
a gravity-mediated type Q-ball. In this case we need to con- 
sider a new type of Q-ball, intermediate between gauge- and 
gravity-mediated Q-balls, whose structure will be determined 
by the region of the potential where the corrections due to the 
messenger scale are still significant. We will study the struc- 
ture of these intermediate Q-balls in a future analysis ll22ll . 

Due to the suppressed A-terms in GMSB models, the en- 
ergy in the condensate fragments is typically larger than the 
energy in the baryon asymmetry by a factor of 10-100. In our 
discussion we have assumed that the condensate fragments of 
charge Q subsequently evolve into Q-balls of the same charge. 
One may question whether this assumption is reasonable. We 
have considered condensate fragments which initially form 
with only a small amount of energy in their baryonic charge 
relative to their mass. This is consistent with the most recent 
numerical simulation of condensate fragmentation, presented 
in [9], where the initial condensate fragments are called "first 
generation Q-balls". Their subsequent evolution is found to 
depend on how ellipitical the condensate trajectory is. It was 
found that if (essentially equivalent to e^^ in |9]) is less 
than about 10, most fragments are quasi-stable and will proba- 
bly slowly evolve into Q-balls. However, if is 100 or more, 
then the initial oscillon-like fragments break up into Q-balls 
of positive and negative charge. Thus our assumption of evo- 
lution to Q-balls is consistent with the results of |9] if ^ 10, 
but in most of our examples is between 10 and 100, where 
the evolution of the fragments is not clear. The non-linear 
dynamics of oscillon-like objects requires very high resolu- 
tion simulations to evaluate their stability [23]. The spatial 
size of the simulations of L9J is sufficiently large to observe 
many condensate fragments, but this necessarily means that 
the resolution of a single fragment is not optimized; a high- 
resolution simulation of just one fragment would be the best 
way to clearly establish the value of r£ for which the frag- 
ments are quasi-stable. The emission of oppositely charged 
Q-ball pairs might be interpreted as a way for the condensate 
fragments to lose excess energy, but this appears to be con- 
tradicted by the observation of quasi-stable fragments with 
rE ~ 10 in the simulation of [9], which suggests complex non- 
linear dynamics favouring spherically symmetric oscillon-like 
states. This is an interesting question for future research. 

Is Q-ball decay in GMSB AD baryogenesis a testable model 
for the origin of gravitino dark matter? The model predicts a 
large messenger scale, which will be reflected in the spectrum 
of SUSY particles |24]. Moreover, the 2 GeV gravitino mass 
predicted by the model is sufficiently large that there should 
be a significant gravity-mediated contribution to the spectrum 
of SUSY particle masses. The model also predicts that decay 
of the thermal relic NLSP density must underproduce grav- 
itino dark matter when mj, is assumed to be 2 GeV, in order 
that the contribution from thermal NLSP decay is less than 
that from Q-ball decay. These features may be testable at 
the LHC; in particular, identification of the 2 GeV gravity- 
mediated contribution to SUSY breaking masses would pro- 
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vide strong evidence^. The GMSB flat direction might also 
produce correlated baryon and dark matter isocurvature per- 
turbations, if the GMSB potential in the angular direction is 
sufficiently flat and |<I>| has the right magnitude during infla- 
tion. 

In summary, we have shown that unstable Q-balls are a 
natural possibility in GMSB AD baryogenesis with a large 
messenger scale. Such unstable Q-balls are essential if AD 

' The phenomenological advantages of an 0(1) GeV gravitino with respect 
to flavour problems and dark matter have been discussed in the context of 



baryogenesis in GMSB is to be consistent with astrophysi- 
cal constraints from neutron star stability. The requirement 
of a large messenger scale is also consistent with the 2 GeV 
gravitino mass required if gravitino dark matter comes from 
Q-ball decay. This points to a striking consistency between 
the requirements for gravitino dark matter from Q-ball decay 
and successful Affleck-Dine baryogenesis in GMSB models. 

'"sweet spot" (25] and "Goldilocks" | 26] supersymmetry. 
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